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In this letter we study how hyperbolic and non hyperbolic regions in the neighborhood of a 
resonant island perform a important role allowing or forbidding stickiness phenomenon around 
islands in conservative systems. The vicinity of the island is composed by non hyperbolic areas 
that almost prevent the trajectory to visit the island edge. For some specihc parameters there are 
tiny channels embedded in the non hyperbolic area that are associated to hyperbolic fixed points 
present in the neighborhood of the islands. Such channels allow the trajectory to be injected in the 
inner portion of the vicinity. When the trajectory crosses the barrier imposed by the non hyperbolic 
regions, it spends a long time to abandon the surrounding of the island, since the barrier also prevents 
the trajectory to scape from the neighborhood of the island. In this scenario the non hyperbolic 
structures are the responsible for the stickiness phenomena, and more than that, the strength of 
the sticky effect. We reveal that those properties of the phase space allow us to manipulate the 
existence of extreme events (and the transport associated to it) responsible for the non equilibrium 
fluctuation of the system. In fact we demonstrate that monitoring very small portions of the phase 
space (namely « 4 x 10“"^ % of it) it is possible to generate a completely diffusive system eliminating 
long time recurrences that result from the stickiness phenomenon. 
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Understanding the transport properties of Hamilto¬ 
nian system is one of the major objective in the statistical 
analysis of dynamical systems [i|, 3. Recent works have 
shown that for a large class of systems, including those 
presenting mixed phase space, the transport can not be 
treated considering just ergodic theory or random phase 
approximation Q. One class of such systems is the low 
dimensional (3/2 or 2 degrees of freedom) Hamiltonian 
system (1/2 degree of freedom corresponds to a periodi¬ 
cal disturbance). Low dimensional Hamiltonian systems, 
commonly present nonuniform phase spaces composed by 
regular (islands) and chaotic regions. The interface be¬ 
tween these regions is far from being a smooth surface 
and the dynamics near the edge between chaotic and 
regular regions is very complex and have been not well 
understood so far. The complexity comes mainly due 
to the presence of stickiness in the boundaries of islands 
The sticky effect forces a trajectory injected into the 
boundary area to stay near the boundary for long periods 
of time. One of the main consequences of this phenomena 
is the existence of power law tails in the Poincare recur¬ 
rence times making the system to present distribution of 
recurrence times displaying algebraically decay for long 
times rather than a exponential decay as expected for a 
normal transport system iilMil- In such a situation 
these systems are characterized as out of equilibrium. 

Some features of the kinetics of Hamiltonian systems 
are important to understand anomalous transport and 
super diffusion. The phase space topology of theses sys¬ 
tems plays a crucial role in the anomalous transport and 
in sticky phenomena 12, 1^. Many problems of science 
such as particle advection in fluids 3, 15 1, transport in 


plasma fusion devices dill , celestial mechanics [3 and 
many others found applications in stickiness occurrences. 

Previous works have studied properties of the bound¬ 
aries between regular and chaotic regions focusin g m ainly 
in the role of the stickiness in the dynamics [y, [l8| and 
also the cantori structures derived from the break of tori 
0. When the system presents stickiness, the dynamics 
of orbits in chaotic sea is observed to be intermittent, 
where after periods of chaotic motion away from the in¬ 
fluence of a sticky island, the system presents periods of 
almost regular motion. However a mechanism based on 
the topology of the island vicinity and how that topology 
affects an injected and subsequently ejected trajectory 
from the sticky area is not completely known. 

Here we study how characteristics of the topology of 
the system, namely hyperbolic and non hyperbolic re¬ 
gions on the neighborhood of a resonant island perform 
an important role in order to establish the presence and 
more than that, the strength of the sticky effect. We show 
that the sticky effect is associated with the presence of in¬ 
jection channels related to the crossing of stable and un¬ 
stable manifolds of hyperbolic fixed points in the vicinity 
of the island, allowing the trajectories to shift between 
sticky and nonsticky areas of the phase space. We show 
that the effectiveness of such channels to capture trajec¬ 
tories to the sticky area is closely related to the degree of 
hyperbolicity of the close surrounding area of hyperbolic 
fixed points located in the neighborhood of an island. Fi¬ 
nally we make use of the presence of the hyperbolic chan¬ 
nels to avoid (control) extreme recurrence events and the 
anomalous transport associated to it, resulting from a 
trajectory injection into the sticky area. 
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Here we characterize a hyperbolic region of the phase 
space S as an ensemble for which the tangent phase 
space splits continuously into stable (SM) and an un¬ 
stable (UM) manifolds. SM and UM are invariant under 
the system dynamics: infinitesimal displacements in the 
stable (unstable) direction suffer exponentially decay as 
time goes forward (backward) [1^. In addition, it is re¬ 
quired that the angles between the stable and unstable 
directions to be uniformly bounded away from zero. In 
this way, in order to quantify the degree of non hyperbol- 
icity related to the phenomena we describe in this letter, 
let us consider an initial condition (po, xg ) and an unit 
vector V, whose temporal evolution is given by 

{Pnt Xn)Vji\\^ ( 1 ) 

where J (p„, Xn) is the Jacobian matrix of the map. For n 
large enough, v is parallel to the Lyapunov vector u(j), x) 
associated to the maximum Lyapunov exponent A„ of the 
map orbit starting by {pg,xg). A backward iteration of 
the same orbit gives us a new vector that is parallel 
to the direction s{p,x), the Lyapunov vector associated 
to the minimum Lyapunov exponent As (2lj| . For regions 
where As < 0 < A„ the vectors u{p,x) and s(p,x) are 
tangent to the UM and SM, respectively, of a point (p, x). 

The (non)hyperbolic degree of a region S can be stud¬ 
ied computing the local angles between the two manifolds 


0(p, x) = cos {\u ■ s|). 


( 2 ) 


for {p,x) G S [^. So, 0{p,x) ^ 0 denotes tangency 
between UM and SM at (p, x). The general method used 
to calculate the 0 angles follows the reference 221. 

Chaotic orbits of two dimensional mappings are often 
non hyperbolic since the SM and UM are tangent in in¬ 
finitely many points. As an illustration of this effect we 
consider a periodically kicked rotor subjected to a har¬ 
monic potential function - the Chirikov-Taylor map Q, 
whose dynamics is two dimensional. The dynamics of a 
periodically kicked rotor can be described in a periodic 
phase space [—7r,7r) x [—7r,7r), whose discrete-time vari¬ 
ables Pn and Xn are respectively the momentum and the 
angular position of the rotor just after the nth kick, with 
the dynamics given by the following equations: 


Pn+i =Pn + K sin(a;„), mod 27r (3) 

Xn+I =Xn+ Pn+1, mod 27r (4) 


where K is related to the kick strength. 

In order to exemplify the dynamics of the Chirikov- 
Taylor map and its sticky phenomena. Fig. [T]-a presents 
a portion of the phase space for a typical trajectory of 
the system for K = 3.0. The denser areas near the is¬ 
land result from the sticky effect due to the time the 
trajectory remains near the edge of the island. In or¬ 
der to characterize the hyperbolicity of the surrounding 
areas of the island Fig. [T]-b displays the local angle be¬ 
tween stable and unstable manifolds (0), Eq. ([21), near 



FIG. 1. (color online) (a) Phase space for the kicked rotor 
map, darker regions around the main island reflects the effect 
of the stickiness, (b) Phase space distribution of angle be¬ 
tween unstable and stable manifold, dark (blue) tones mean 
strong non hyperbolicity. (c) An example of the edge detec¬ 
tion algorithm used here to compute the vicinity of the island. 
In the inset we display a magnification showing details of the 
edge, (d) Probability distribution function of the angle be¬ 
tween stable and unstable manifolds, p{0)- 


the resonant island. It is clear that the major part of the 
vicinities of the island is composed by strong nonhyper- 
bolic areas, represented by dark blue areas in Fig. [T]-b 
(tangencies between UM and SM). A small fraction of 
the vicinity is observed in red - yellow - green tones and 
corresponds to angles greater than 30° and responsible 
for the weak hyperbolic part of the vicinities. The role 
of both areas will be clear later on in the text. In Fig 
mh we also define the angle (f) = arctan(p/a;), defined in 
the interval [—7r/2,-|-7r/2], so each point in the vicinity 
of the island can be identified by a single number. The 
9-1-9 (symmetric) hyperbolic fixed points present in the 
neighborhood of the main island are identified as black 
bullets in Fig. |T]-b. In order to make clear what we call 
the vicinity of the island. Fig [T]-c displays the result of 
our algorithm for edge island detection, details magnified 
in the inset. Finally in Fig[T]-d we plot the probability 
distribution function of the angles between SM and UM 
p[9), so p{0)d0 represents the probability to find a angle 
between 9 and 9 + d9 m the phase space ensemble dis¬ 
played in Fig. |T]-b. The large plateau for small 9 angles 
reflects the strong nonhyperbolic character of the region. 

The topological properties of the phase space in the 
vicinities of an island play an important role in the sticky 
mechanism. In order to explore the relation between 
topological characteristics of the phase space and the way 
trajectories visit an island vicinity and stick to it, we de¬ 
fine the probability density function for trajec- 
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FIG. 2. (a) Probability distribution function of income angle 
in the vicinity of a island of the Chirikov-Taylor map. (b) 
Average time spent in the sticky area as a function of the 
injected angle, (c) Magnification of panel (a) and (b) near the 
maxima of {(p) 
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tories injected into the vicinities of the island consider¬ 
ing the vicinity computed by our algorithm (Fig. [T]-c). 
F^^\(p)d(j) is the probability that a typical chaotic tra¬ 
jectory to visit the vicinity of the island through a angle 
between (j) and (/) -|- d(j). In Fig. [2] we plot (</>), panel 
(a) as well as the average time the trajectory remains 
in the vicinity of the island when injected by a particu¬ 
lar angle, panel (b). In panel (c) we magnify the gray 
region of panel (a) and (b). Observe that although the 
major part of the trajectories are injected by just few an¬ 
gles inside the red-yellow-green tones regions around the 
island in Fig. [T]-a (weak hiperbolic regions ), these spe¬ 
cific trajectories spend, in average, a short time mapping 
the sticky area. Trajectories are easily injected into the 
sticky areas by weak hyperbolic areas surrounding the is¬ 
land but almost all of them are also easily ejected from it. 
Those trajectories do not contribute for the phenomena 
of stickiness and do not make any substantial changes in 
the Poincare recurrence time for the dynamics. 

In order to distinguish sticky trajectories from those 
that just reach the island edge and leave it quickly, we 
compute the probability density function (</)), of 

trajectories injected into the sticky area by a specific an¬ 
gle considering that once a trajectory reaches the vicinity 
of the island, it remains mapping the same set of points 
as computed by our algorithm of island edges detection 
for at least 100 iterations. We identify such trajectories 
as sticky ones. The result is plotted in Fig. [3]-a. Almost 
all sticky trajectories are injected in very specific inter¬ 
vals of angles. Each angle intervals are directly related to 
the angular location of the chain of periodic points (set as 
black bullets in Fig. [T]-c). It is possible to conclude that 
the trajectories are injected into the sticky area when 
they are tangent to the stable manifold of the 18*^ or¬ 
der fixed points chain located in the vicinity of the main 
island. All trajectories that are not tangent enough to 


(t> 

FIG. 3. (a) Probability distribution function of income angle 
into the vicinity of the main island of the Ghirikov map. (b) 
Average spent time in the sticky area as a function of the in¬ 
jected angle, (c) Probability distribution function of outcome 
angle from the vicinity of the main island of the Chirikov map. 
(d) Average spent time in the sticky area as a function of the 
ejected angle, (e) Probability distribution of injected angle, 
subjected to the condition that the time of stickiness to be 
greater than 1000 iterations 


the stable manifold of the hyperbolic fixed point do not 
cross the tiny hyperbolic channel produced by the cross¬ 
ing between stable and unstable manifolds of the fixed 
point and can not be injected into the stick area. 

Considering those trajectories injected into the stick¬ 
iness remaining mapping the edge for at least 100 iter¬ 
ations, we compute in Fig. [5]-b the average time they 
spend near the island (sticky trajectories) as a function 
of the injected angle. In Fig. [3]-c we graph the prob¬ 
ability distribution function of sticky trajectories as a 
function of the ejected angle Fq™((^). It is clear the al¬ 
most discrete nature of the distribution. All ejected tra¬ 
jectories follow the unstable manifold of the hyperbolic 
fixed points moving along a narrow channel departing 
from the fixed point. We plot in Fig. [3]-d the average 
time the trajectories stay in the sticky region (at least 
100 iterations inside the sticky area) as a function of the 
outcome angle. From Fig. [5]-d we can conclude that 
sticky trajectories are ejected only by few angle intervals 
4>out = ?f>(max(roKt)). Therefore, we are able to calcu¬ 
late the probability J-{(j))d(j) that a given trajectory will 
enter the sticky region considering only trajectories that 
leave these region through the angle (font- The result is 
plotted in Fig. |31-e. The great similarity between T and 
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suggests that the or previous conclusions are con¬ 
sistent. Therefore, we can argne that the local maximum 
of represent the sticky angles, i. e., the angles that 

once a trajectory is injected from one of them, there is a 
great probability that this trajectory turn to be stuck to 
the island. These maxima correspond to the same region 
where are located hyperbolic points, confirming the hy¬ 
potheses that these points provide channel for a typical 
trajectory to enter in the sticky region. Figs. [31-a and 
[3]-e, clearly show us that both figures are almost iden¬ 
tical, a strong suggestion that all trajectories leave the 
sticky regions by the hyperbolic channels departing from 
the hyperbolic fixed points. 

The presence of stickiness around an island is predicted 
by some theory [l^ but an analysis of Fig. [3]-a shows 
that the effectiveness of a trajectory injection or ejection 
by a particular hyperbolic channel is not the same for all 
fixed points, as can be observed by the different ampli¬ 
tude of maxima of in Fig. [3]-a. In order to make 

clear the role of the hyperbolicity of the close vicinity of 
the hyperbolic fixed points in the injection and ejection 
phenomena of sticky trajectories, firstly, we present as an 
example, in Fig. |3]-a the degree of hyperbolicity of one 
of the 18 fixed points presents around the main island 
of the Chirilov-Taylor map. As observed in Fig. 0]-b the 
probability density function p{0) presents just a sharp 
maximum due to the almost unique angle between sta¬ 
ble and unstable manifold computed in the close vicinity 
of the fixed point. All other fixed points present simi¬ 
lar sharp peaks in the probability density function p(0), 
nevertheless each hyperbolic point has its own angle for 
the maximum of p{0) characterizing its own degree of 
hyperbolicity. Secondly, to demonstrate the relation be¬ 
tween the degree of hyperbolicity of the close vicinity of 
the fixed points, and the effectiveness of the hyperbolic 
channels related to each fixed point to capture sticky tra¬ 
jectories, we graph in Fig. [5] the maxima of the function 
as a function of the degree of hyperbolicity mea¬ 
sured by the angle for which the function p{0) presents 
a maximum 0(pmax)- The red line is a power law fit¬ 
ting that serves us as eye guide. The result presented in 
Fig. [S] clearly shows that the effectiveness of the chan¬ 
nels is a function of the degree of hyperbolicity of the 
close region of the fixed points. Small values of d(pmax) 
is related to the fact that just a very small portion of 
the surrounding area of the fixed point is occupied by 
the injection/ejection hyperbolic channel. As a result 
the function presents a relatively small maxi¬ 

mum, meaning that just a small fraction of trajectories 
can cross the channel in an injection or ejection process 
from the sticky area. 

Such properties of the phase space allow us to manip¬ 
ulate the non equilibrium fluctuation of the system. To 
show that it is possible to control the non equilibrium 
fluctuations that arise due to the presence of stickiness. 



FIG. 4. (color online) (a) Degree of hyperbolicity of the vicin¬ 
ity (e = 0.0025) of one of the eighteen fixed points present 
around the main island of the map U) (b) p{9) character¬ 
ized by just on maximum when computed near a fixed point 
around the main island 



FIG. 5. (color online) The effectiveness of the injection chan¬ 
nels measured using the maxF)^° as a function of the angle 
between UM and SM. 


we track the position of the trajectory and once it maps a 
small circle of radius 0.003 centered in one of the 18 fixed 
points we perturb the trajectory, so a possible crossing 
of the channel and consequent stick of the trajectory is 
avoided. Numerically, we perform a restart of the tra¬ 
jectory outside the injection channel. Results for the 
Poincare recurrence time for the system with and with¬ 
out the control mechanism for two values oi K, K = 3.0 
and K = 3.565 (a large stickiness case) are plotted in 
Fig. [6l Black bullets and green up-triangle display the 
time distributions for Poincare recurrence without any 
control mechanism. As can be observed for large recur¬ 
rence time, a strong fluctuation of the exponential law 
is observed. In fact for large recurrence time the dis¬ 
tribution has a power law decay as a result of the sticky 
phenomena in the recurrence time. The time distribution 
of the Poincare recurrence time for the system subjected 



FIG. 6. (color online) Recurrence time for the Chirikov-Taylor 
map with and without control 
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to our control mechanism is displayed red squares and 
blue down-triangle. For this case, almost all fluctuation 
for long recurrence time is absent, corroborating to the 
idea that all nonequilibrium fluctuation in the system is 
now absent since the stickiness is avoided. Additionally 
observe that the exponential rate for both K values is 
the same supporting the idea that the behavior of the 
system is now completely diffusive independently of the 
K value. All results here are presented for two values of 
K but similar results are obtained for other values of the 
nonlinear parameter. 

In conclusion, we describe in this letter a mechanism to 
suppress the effect of stickiness based on the knowledge 
of the nonhyperbolic structure on the edge of a island 
of a Hamiltonian 2 degree of freedom system. We show 
that the effectiveness of an island edge to stick trajecto¬ 
ries is directly related to the degree of hyperbolicity of 
small areas surrounding fixed points around the island. 
We show that monitoring those areas of the phase space, 
it is possible to generate a complete diffusive processes 
without (almost) any influence of the large recurrence 
time due to the stickiness phenomena. Since very large 
sticky times make the system to present extreme events 
in the dynamics, it possible to affirm that, once under 
control, we can turn the out-of equilibrium system into 
a in equilibrium one. 
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